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Abstract 

Ultrahigh energy protons in magnetic fields produce pions and thus 
lose energy. The mean free path of such a process is worked out for 
Gaussian random fields. Two cases are considered: an isotropic and a 
cylindrically symmetric distribution. The energy loss is proportional 
to E^(B 2 ); it becomes significant for protons of energies ^ 10 19 eV 
and magnetic fields B ■> 10 9 Gauss. For energies and magnetic fields 
of this magnitude, a proton injected into the magnetic field loses a 
substantial fraction of its initial energy due to pion production. 

PACS: 13.85.Tp, 98.70.Sa, 98.70Rz 



1 Introduction 

It is well known that the propagation of ultra high energy (UHE) protons in 
the Universe is limited by the Greisen-Zatsepin-Kuzmin (GZK) mechanism, 
|], ||. The propagating protons undergo inelastic scattering on the photons 
of the cosmic microwave background radiation (CMBR) and produce pions. 
On the average, the initial energy is shared (roughly) equally by the nucleon 
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and pion in the final state. The energy at which the GZK mechanism be- 
comes important can be crudely estimated by saturating the inelastic cross 
section by the A resonance. Given the fact that the average energy of the 
CMBR photons is around 3 x 10~ 4 eV, one gets that pion production becomes 
significant around proton energies of the order of 10 19 eV. (In fact, a similar 
simple estimate was used in Greisen's original paper.) 

Here we discuss a mechanism for energy loss of UHE protons, hitherto ap- 
parently ignored, viz. by inelastic scattering on virtual photons. This mech- 
anism plays no significant role in limiting the propagation of UHE protons in 
intergalactic space. However, it becomes significant when the propagation is 
considered in an environment where strong magnetic fields exist, e.g. in jets 
emerging from gamma ray bursters (GRB) or jets in active galactic nuclei 
(AGN). 

The "average" energetics of pion production on virtual photons (typically, 
on an external magnetic field) is very different from that of the GZK mech- 
anism. In fact, if the size of the external magnetic field is characterized by 
a length L, then the typical momentum of the virtual photon is of the order 
of 1/L. Consequently, the average invariant center of mass (CMS) energy 
available for pion production is of the order, 

s ^m 2 + 2E/L, 

where m stands for the nucleon mass and E is the energy of the incident 
proton in the rest frame of the local universe. Obviously, this energy is much 
less than the analogous quantity in the GZK process for any macroscopic 
L. Nevertheless, for a sufficiently strong and spatially confined magnetic 
field, one obtains an appreciable production rate, due to the fact that the 
Fourier spectrum of a confined field is rather slowly decreasing with the wave 
number (typically, as a power). As a consequence, Fourier components with 
|k| 3> 1/L can play a significant role. 

This paper is organized as follows. In the next section, we obtain a 
general expression of the cross section for the process p + B — > X, where B 
stands for an external magnetic field. We also discuss the approximations 
one can make in order to simplify the calculation. The subsequent section, ^ 
contains an evaluation of the interaction rate for random magnetic fields: we 
believe that this serves as a first model for energy loss in the chaotic fields 
present in typical astrophysical environments. Two situations are considered 
in detail: an isotropic and a cylindrically symmetric probability distribution 
of the random field. 
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While an isotropic environment largely serves to illustrate the physical 
features of the process in a simple context, it is also potentially applicable to 
a situation in which the size of the magnetic field is substantially larger than 
the interaction mfp. The calculation of the interaction rate for a cylindrically 
symmetric environment is relevant for jets, for instance, emerging from an 
AGN or GRB. The results are discussed in sec. |j. 



2 General expression of the cross section in 
an external magnetic field. 

The calculation described here is an elementary application of the optical 
theorem. The amplitude of a proton interacting with an external field and 
producing a final state \X) is given by 

T(p + B^X) = Jd 4 x A lt (x)(X\j"(x)\p) (1) 

Here, A^(x) stands for the vector potential of the external field and j^x) is 
the density of the electromagnetic current. Squaring (|l|) and summing over 
the final states \X) , one expresses the cross section in terms of the current 
correlation function. This is textbook material, see for instance 0. We 
also assume that that the external magnetic field is static. This assumption 
simplifies the calculation. From a physical point of view, it is justifiable even 
if astrophysical objects of bulk Lorentz factors of the order of a few hundred 
are considered: the protons we are interested in have Lorentz factors which 
are ten or eleven orders of magnitude larger. 

On writing for the Fourier transform of the vector potential 

MO) = S (?o) % (q) (2) 
and using a gauge in which A = 0, one gets: 

o = ^^\d\ a*(q)a k (q)W ik (3) 

In equation (^) m and E stand for the mass and energy of the incident proton, 
respectively and Wik is the spatial part of the standard polarization tensor: 

u/ Fl ( j_ QtfA , F 2 / (pq)\ ( (pq)\ m 
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The notation is standard, p and q are the four momenta of the incident proton 
and virtual photon, respectively, v = (pq) / m. 

A further simplification is possible due to the fact that the protons we are 
interested in are extreme relativistic and the average value of the momen- 
tum of the virtual photon is of the order of 1/L. In order to motivate this 
simplification, we Lorentz transform to the rest frame of the proton. In that 
frame the components of the four momentum q and the field quantities are 
distinguished by a prime. Components perpendicular to the direction of mo- 
tion are denoted by capital letters; longitudinal components by a subscript 
/. Since we have v ~ 1, the transformation formulae are: 

%~\ e *-V{y)<lu q'i ~ ^exp(y)q h q A ~ q A . (5) 

B l = Bi, E l — Ei — 0, 
b 'a ~ 2 ex P(y) B A, E' A ~ - exp(y)e AB B B . (6) 

In the last two equations, y stands for the rapidity. 

As a consequence, apart from corrections of 0(exp(— y)), 

q 2 ~ 0, B ■ E ~ 0, B 2 — E 2 ~ 0. (7) 

In a reference frame comoving with the proton, the magnetic field appears 
as a stream of (almost real) photons: consequently, the contribution of the 
structure function F 2 to the cross section is negligibly small. 

One can then express the cross section on the external magnetic field in 
terms of the photoproduction cross section, <r 7 , viz. 

a ~ ^ J d 3 q fli(q)aj(q)* P ' q <7 7 (^q 2 - qiq 3 ) (8) 

One readily recognizes that eq. @ is equivalent to a Weizsacker- Williams 
approximation to the cross section. Because of the presence of a transverse 
projector, that expression is a manifestly gauge invariant one. It is worth 
noticing that in the Weizsacker- Williams approximation the expression of 
the cross section is independent of the mass of the projectile. Hence the 
same expression can be used to describe e.g. photon induced reactions in a 
magnetic field. 
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Finally, one considers the evaluation of a T Due to the fact that the 
photoabsorption cross section is to be evaluated near the pion production 
threshold, to a good approximation one can saturate it by the contribution of 
the A resonance. A narrow resonance approximation is sufficiently accurate. 
Hence we put 

ct 7 w a 5 (s - rn^) , (9) 

where the dimensionless quantity, a is the integral of the pion photoproduc- 
tion cross section across the resonance, 




a(s)ds. 



Using a standard invariant Breit-Wigner fit and the data available, Q one 
gets cr ~ 0.3. 



3 Random magnetic fields 

We model the chaotic magnetic fields present in the astrophysical environ- 
ments of interest by means of a Gaussian random field of zero mean. The 
central object in the theory of random fields is the generating functional 
of the correlation functions. In the case of a Gaussian field, only the sec- 
ond cumulant is different from zero. We write the generating functional as 
follows IHl. 



Z[j] = J Vsl exp -S + i J d 3 kj r (k) a r (k) 



(10) 



In eq. (|10J) j stands for an external source, a is the Fourier transform of the 
vector potential, cf . eq. (|2|). The functional S is a generalized entropy; for 
a Gaussian field it is a (gauge invariant) quadratic functional of a. We write 
S as follows. 



S = Jd 3 k a*(k) aj (k) - -U. J 4 (l + L 2 k r k s u rs ) . (II! 



In eq. (0), L stands for the root mean square correlation length (the av- 
erage is taken over directions). The tensor u TS characterizes the directional 
distribution of the probability density. For a general, arbitrarily anisotropic 
distribution, u TS has 6 independent components: e.g. the three, mutually 
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orthogonal, principal correlations and the three angles describing the orien- 
tation of the principal correlations. In what follows, however, we consider 
environmernts of high symmetry; consequently, fewer parameters are suffi- 
cient. The factor (1 + L 2 k r u rs k s ) ensures an exponential decrease of the 
correlation function with distance, cf . ref. |J. 

In considering particle production in a random field, one has to replace 
factors such as a*a,j in eq. (|3|) and subsequent ones by their expectation values 
in the ensemble defined by eqs. (fLOD and (|TTf) . 

We now consider two special cases of the ensembles in order to calculate 
particle production cross sections. 



3.1 Isotropic ensemble. 

This ensemble is characterized by the tensor Uij in eq. flTT| ) being the unit 
tensor, = 5ij. One finds: 

(a 4 ( q) a j(q r)=5 3 ( q - q ) (^-^)^f^(l + ^qT 2 (12) 

In the expression of the cross section, however, one finds a factor (ai( q )a*( q )), 
which is infinite, see the last equation. This is due to the fact that we idealized 
a region of non vanishing magnetic field by one of infinite extent albeit of 
an exponentially decreasing correlation function. In order to correct for the 
inconsistency caused by this idealization, we replace the delta function of 
vanishing argument by a quantity proportional to the volume, viz. 

5 3 (0) — > -V. 
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(This is a consistent procedure provided the density of levels can be approx- 
imated by the Rayleigh - Jeans formula, as done here. In the problem under 
consideration, the conditions for the validity of that approximation are sat- 
isfied.) We take V to be the correlation volume; thus, in the isotropic case, 
V = 47rL 3 /3; clearly, different geometries give rise to different expressions 
of the correlation volume. The important fact is, however that the incident 
flux is oc 1/V. Thus the reciprocal mfp is independent of the choice of the 
volume. 

With this and using eq. @ the cross section can be evaluated in terms 
of elementary functions. Quoting directly the inverse of the absorption mfp 
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which is the relevant quantity for the applications, one finds: 

J_ = L (r2) vq 
A„ 



7T 



m. 



r/H- 



(13) 



Here / is a function of the dimensionless variable, w = {m\ — m?) L/2E. Its 
explicit form is: 

'1 + w T 



f(w) = w 
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+ 



1 + w 2 



(14) 



However, it was pointed out in the Introduction that 1/L is a small momen- 
tum. As a consequence, we only need eq. (H) for large values of w. In that 
case, eq. ([14]) simplifies to: 



2w 2 



Hence, the expression of A a becomes: 



{w > 1). 



E 2 



A 



Ln 



mr 



(15) 



3.2 Cylindrically symmetric ensemble 

This geometry is a more realistic one. In particular, an astrophysical jet as 
emerging, for instance, from a GRB or an AGN can be approximated by 
a cylindrical geometry at the early stages of expansion. (At early stages, 
the lateral expansion is negligibly small compared to the longitudinal one.) 
Approximating the jet by one of cylindrical geometry means that the lateral 
expansion is neglected altogether. It has been known for a long time that 
this is an acceptable approximation in the initial stages of expansion of a 
relativistic fluid 

In this case, the tensor Uij in eq. (|TT|) effectively depends on one parameter 
only. It is convenient to introduce the longitudinal and transverse correlation 
lengths with respect to the axis of the cylinder and an anisotropy parameter, 
a, such that 

L 2 T = aL 2 , L 2 L = L 2 {l-a). 
In practice, «<1, say a ~ 0.1 or soQ Using this parametrization, we have: 

(l + L 2 qi u ijqj ) = (l + (aq T 2 + (1 - a)qfj) (16) 
1 E. Vishniac, private communication. 
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In eq. (IB), qx and qi stand for the momentum components perpendicular 



and parallel to the axis of the cylinder, respectively. 

In the case of such a geometry, the integral occurring in eq. (|8[) cannot 
be calculated in a closed form. In essence, this is due to the fact that the 
expression of the absorption cross section now contains two directions: that 
of the incident proton and the axis of the cylinder. However, instead of 
resorting to a numerical evaluation, we observe that the variable w is large 
and the absorption can only be significant if the angle between the incident 
proton and the axis of the cylinder is not too large: efficient absorption 
requires a coherent magnetic field. 

These simplifications allow a calculation of the absorption mfp in a closed 
form. A somewhat tedious, but elementary calculation leads to the result: 



1 4 ff0 E 2 (B 2 ) 

7 T b l B J^ ~ 



The factor F(Q) is given by the expression: 

F{Q) « (cos 9) 3 In ^ (cOS a 9)2 - 1^ (18) 

In eq. flIB|) stands for the angle between the incident proton and the axis 
of the cylinder. Obviously, this expression holds only if the angle is small. 
From the physical point of view, however, this is not a serious limitation: due 
to the presence of the factor oc (cos Oy, the mean free path becomes very 
large unless the angle of incidence with respect to the axis of the cylinder is 
small. 



4 Discussion 

Our approach has the advantage that it does not depend on the details of 
the production process, since it is based on the use of the optical theorem. 
However, its limitation is that the absorption cross section is obtained to 
lowest order in the fine structure constant. This poses no problem as long as 
J (B 2 ) < m 2 / e = -B cr jt, where m is the mass of a charged particle involved in 
the process. For electrons and light quarks (u,d), the value of B CT ^ is around 
10 14 Gauss. In magnetic fields of this order of magnitude, radiative corrections 
and pair production become important. To our knowledge, no results are 



8 



available for such field strengths. Existing calculations, such as Erber's, J7J 
assume a homogeneous magnetic field. Calculations of this type can be used 
to estimate energy losses as long as the magnetic fields are approximately 
homogeneous on the scale of the Larmor radius of the propagating charged 
particle. For realistic circumstances, however, this is hardly the case. Thus, 
the question about the energy loss of charged particles in astrophysically 
important magnetic field approaching the critical value of the field, is still an 
open one. Our formulae, however, are expected to give at least a qualitative 
insight into the question of absorption even for near-critical fields. 

Our results show that the circumstances neeeded for the applicability of 
eq. (p~5| ) are hardly met: one needs magnetic fields with a coherence length 
substantially in excess of the Larmor radius at high energies. Nevertheless, 
that equation is an instructive one: due to its simplicity, the general features 
of the absorption cross section is easily understood. 

From the physical point of view, eq. (|17D is more interesting. In order to 
assess the importance of the process discussed it is worth converting eq. (|17| ) 
into a form permitting numerical estimates. The value of <To has been quoted 
before; the rest of the numbers is also taken from ref. 0. One obtains: 

X a L [U) llCPeW (lO^Gauss) 2 1 J 

(We used the usual conversion factor between the natural and conventional 
units of the magnetic field, viz. B/l(MeV) 2 = 1.9 x 10~ 14 .B/lGauss.) 

We find that a paraxially propagating proton of energy ~ 10 20 eV travers- 
ing a - relatively modest - magnetic field of 10 9 Gauss has an absorption mfp 
about (l/5) th the size of the magnetic field. 

In a collision at the relevant energies, on the average the nucleon and 
the produced pion in the final state share the incident energy equally. As 
a consequence, using the continuous energy loss approximation, the energy 
loss per unit path length is: 

In the last equation, e stands for the fractional energy loss of a nucleon (in the 
laboratory system) due to pion production. Assuming as we do throughout 
this paper that the pion production cross section is dominated by the A 
resonance, one gets, 

ml 

e — 

2m&m 
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Due to the fact that 1/A oc E 2 , the energy loss per unit path length grows 
as E 3 . 

In fact, by inserting the expression for the mfp given by eq. (|T9| ) into 
eq. (|20D, the equation for the energy loss is readily integrated. We exhibit the 
result for the energy loss of paraxial protons (F(Q) ~ 1) over one correlation 
length, L. We get: 



E(x) 
W) 



1 + 0.1 



E \ 2 / (B 2 ) \ j 



10 20 GeV/ V(10 9 G auss / ^ 



-1/2 

(21) 



We used the mass values listed in ref. || in order to arrive at eq. (|2T|). 

We conclude that the mechanism described in this paper appears to be a 
major obstacle to accelerating protons up to energies of the order of 10 19 eV 
or more by a conventional Fermi acceleration mechanism. One notices for 



instance that a proton of E — 10 20 eV injected into a field of y (B 2 ) = 
10 10 Gauss loses about 70% of its initial energy over a correlation length. 

This adds to the puzzle of the highest energy cosmic rays: it is known 
that particles of energy about 10 20 eV arrive to the Earth and they give rise 
to extensive air showers. At the same time, it appears to be increasingly 
difficult to find an efficient mechanism for producing them at the usually 
suspected sites, for instance in active galactic nuclei or gamma ray bursters. 
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